Classical analysis of correlated multiple ionization in strong fields 
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We discuss the final stages of the simultaneous ionization of two or more electrons due to a 
strong laser pulse. An analysis of the classical dynamics suggests that the dominant pathway for 
non-sequential escape has the electrons escaping in a symmetric arrangement. Classical trajectory 
models within and near to this symmetry subspace support the theoretical considerations and give 
final momentum distributions in close agreement with experiments. 
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I. INTRODUCTION 

Classical models for atomic processes can provide use- 
ful insights in situations where quantum effects are not 
prominent, as for instance in the dynamics of highly ex- 
cited states or in multiphoton absorption processes. Mi- 
crowave ionization of hydrogen and other atoms with a 
single valence electron localized wave packet dy- 

namics or electron scattering off atoms come to mind 
They provide a natural starting point for semiclas- 
sical investigations that include, at least approximately, 
quantum interference effects. In particular the ioniza- 
tion of Rydberg atoms in linear and elliptically polarized 
microwave fields has received considerable attention and 
similarities and differences between classical and quantal 
behaviour have been sorted out in great detail In 
most cases only a single electron is considered and suffi- 
cient to interpret the observations. 

Interactions between electrons seem to play an impor- 
tant role in multiphoton multiple ionization in strong 
laser fields. Experiments show that the yield of multi- 
ply charged ions is much higher than can be expected on 
the basis of an independent electron model |^,||. More 
recently, it has been noted that the electrons can escape 
non-sequentially and that they are correlated in their fi- 
nal state ]7|~p^. This correlation in the final state came 
as a big surprise and it is our main objective here to 
discuss a classical model for it. 



II. THE MODEL 

The full process of multiphoton multiple ionization is 
quite complicated and involves many steps. A plausible 
model relevant for the field intensities of the experiments 
is the rescattering mechanism pd]-[l^. Before the pulse 
arrives, the atom is in its ground state. Then one electron 
escapes from the atom, most likely by tunnelling through 
the Stark barrier. This electron is then accelerated by the 
field and can be reflected back towards the atom. During 
this impact energy is transfered to other electrons, per- 
haps lifting them above the ionization threshold or bring- 
ing them close enough so that tunnelling is again possible. 
If not enough energy is provided at this stage, perhaps 
another rescattering process can follow until eventually 



multiple ionization takes place or the pulse disappears. 
However, before the escape to multiple ionization all ex- 
cited electrons pass close to the nucleus where they inter- 
act strongly with the each other and with the Coulomb 
attraction. During this phase their (classical) motion is 
very fast compared to the field oscillations and an adi- 
abatic analysis, keeping the field fixed, can be applied. 
Morevoer, because of the strong interaction all memory 
of the previous motion is lost, so that the initial state 
for the multiple ionization event is a statistical distri- 
bution of electrons close to the nucleus. Our discussion 
starts once this intermediate cloud of excited electrons 
has formed. We do not consider the process by which 
it has been generated; for instance, one might imagine 
exposing ions to both an electric field and an electron 
beam. 

The arguments that follow focus on two electron escape 
but can easily be extended to discuss multiple ionization, 
as indicated below. The Hamiltonian then has as usual 
the kinetic energy of the electrons, their mutual repul- 
sion, the attraction to the core and the potential due to 
the oscillating electric field. In many experiments the 
recoil momentum of the ion is measured, and given the 
extremely small momenta of the photons it is possible to 
calculate it as the sum of the momenta of the electrons. 

Initially, there is no field and the atom is in its ground 
state. In the final state, after the pulse is turned off, 
both electrons are free and have positive total energy. 
Not all the energy difference has to be provided by the 
impacting electron since there can be additional accel- 
eration by the field after the electrons escape from the 
core region. However, within the adiabatic assumption 
motivated above, the energy content of the intermediate 
electron cloud has to be high enough to let both elec- 
trons escape from the immediate vicinity of the nucleus. 
Without field this implies positive energy, but if the field 
is on and non-zero, a Stark saddle forms some distance 
away from the nucleus and the electrons can escape over 
it. The rapid acceleration downfield will then pull the 
electrons out and feed in the energy needed to remain 
asympotically free once the pulse is gone. 

The Stark saddle that forms in the field provides a fo- 
cus and a bottleneck for the electrons which they have 
to cross in order to leave the atom. All electrons see 
the same saddle and would like to cross it, but if they 
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try to do simultaneously, as suggested by the observed 
electron correlations, their mutual repulsion gets in their 
way. Suppose that one electron is slightly ahead of the 
other when running up the hill towards the Stark sad- 
dle: the one that is ahead has the advantage that the 
repulsion with the companion pushes it uphill, whereas 
the one behind not only has to fight the attraction to the 
nucleus but also the repulsion from the one ahead. Their 
interaction is perfectly balanced if they cross the saddle 
side by side, with reflection symmetry with respect to 
the field axis. The previous considerations suggest that 
deviations from this symmetric configuration are ampli- 
fied and cannot lead to simultaneous ionization. The 
arguments used here are similar to the ones advanced by 
Wannier in his analysis of double ionization upon elec- 
tron impact [|9|j20|. 

Therefore, we propose that near the threshold for 
double ionization the dominant path leading to non- 
sequential double ionization has both electrons escape 
symmetrically with respect to the field axis. If more 
than two electrons are ionized simultaneously the natu- 
ral extension is that they form a regular n-gon in a plane 
perpendicular to the field axis. 



III. SYMMETRIC DOUBLE IONIZATION 

With the field pointing along the a;-axis and two elec- 
trons confined to the plane 2 = their coordinates in 
the symmetric subspace are {x, y, 0) and (x, — y, 0) in po- 
sition and {px,Py,0) and (p^, —py,0) in momenta. The 
classical Hamilton function for this geometry then is (in 
atomic units) 



H{px,Py,x,y,t) ^pI+p'1 + V{x,y,t) 



(1) 



with potential energy 
4 



V{x,y,t) = - 



+ — + 2Fxf{t) cosfwt + . 



and the pulse shape 

fit) ^ sin" (TTt/Td) 



(2) 



(3) 



where the duration of the pulse is taken to be four field 
cycles, Td — Stt/uj. The rescattering of the electrons leads 
to a highly excited complex of total energy E which ev- 
ery now and then is close to the symmetric configuration 
described by the Hamiltonian (1). Any configuration on 
this energy shell (for some fixed time t) as well as any 
phase (j) of the field is equally likely, and the experimen- 
tal observations are averages over initial conditions and 
phases. 

Fig. |l| shows equipotential lines for the potential (||) at 
a maximum of the field for F = 0.137 a.u., corresponding 
to an intensity of 6.6 • 10^"* W/cm^. The saddle is located 



along the line x = rg cos 6 and y = rg sin 6* with 6 = n/Q 
or 57r/6 and at a distance = •\/3/|F/(t) cos{ujt + 0)|. 
For the above mentioned field the saddle has an energy 
of K = -1.69 a.u.. 




FIG. 1. Adiabatic potential V{x,y,t) for fixed time t in 
the symmetric subspace. The saddle moves along the dashed 
line. 
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FIG. 2. Trajectories in the symmetric subspace. One 
frame shows the time evolution of the radial distance r{t) 
(continuous line) together with the instantaneous position of 
the barrier (dashed line) and the other shows the energy. 

Sample trajectories within the symmetric configura- 
tion are shown in Fig. ||. It is evident that they cross the 
saddle during a maximum of the field and that once on 
the other side the energy increases rapidly. This accel- 
eration is accompanigned by a rapid separation from the 
nucleus, so that during field reversals the electrons will 
not return to the nuclues and will remain essentially free. 

A quantity of direct experimental interest is the dis- 
tribution of momenta of the ion, estimated as pion ~ 
^(Pi + P2) • The many reafization of the exper- 
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iment and the unknown details of the initial conditions 
can be modelled by averaging over all initial conditions of 
prescribed energy and all phases of the field. The results 
for peak field strength F = 0.137 a.u. are shown in Fig. ||. 
For initial energy E = —0.58 a.u. the final distribution of 
momenta clearly shows the double hump structure also 
seen in experiment. The perpendicular momentum of a 
single electron show a clear suppression near zero momen- 
tum and a long tail. The suppression for small momenta 
is a clear sign of electron repulsion. We are not aware of 
expeperimental data on this distribution. 
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FIG. 3. Final momentum distributions for F — 0.137 a.u. 
and initial energy E = —0.58 a.u.: (a) ion momentum paral- 
lel to the field axis and (b) perpendicular momentum of one 
electron. 



IV. THE SADDLE AND THREE-DIMENSIONAL 
MOTIONS 



ries started slightly outside the symmetry plane (Fig. |^). 
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Within the symmetric subspace mentioned before the 
position of the saddle separating trapped motion from 
ionized motion is clear. And as in many models of chem- 
ical reactions it has one unstable direction that defines 
the reaction coordinate and a stable motion perpendic- 
ular to the reaction coordinate. However, in the space 
of six degrees of freedom of the full 3-d two electron mo- 
tion and in the adiabatic approximation for the field the 
stability analysis of the saddle reveals two unstable di- 
rections and four stable ones. The stable directions are 
of minor importance: if excited they persist in the neigh- 
borhood of the saddle as some uncoupled finite amplitude 
motions. The second unstable direction besides the reac- 
tion coordinate is responsible for the importance of the 
symmetric subspace. Motion leading away from this sym- 
metric subspace will typically have one electron escaping 
and the other returning to the nucleus. This corresponds 
to single ionization. The electron returned to the nucleus 
may have enough energy to ionize in the next step or may 
gain additional energy from the field to ionize later. Ei- 
ther way, the electrons do not escape symmetrically and 
simultaneously, so that there are no correlations between 
the two outgoing electrons and the ionization is sequen- 
tial. 

Without going into the technical details of this analy- 
sis, we can illustrate some of these features with trajecto- 
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FIG. 4. Trajectories of electrons outside the symmetric 
subspace for E — —0.58 a.u.. Initial positions are close to 
the saddle and marked by heavy dots; the electrons are dis- 
tinguished by dotted and continuous tracks. Frame (a) shows 
for reference the ionization in the symmetric subspace. Frame 
(b) shows a case where outside this symmetric subspace one 
electron escapes and the other falls back to the ion. Frame (c) 
shows an example of sequential ionization of both electrons 
in opposite directions. 

Fig ^a shows initial conditions on the saddle and 
symmetrically escaping electrons. For some deviation 
from symmetry, one electron escapes, the other remains 
trapped to the nucleus (Fig. ^). It is possible, how- 
ever, that the second electron picks up enough energy to 
ionize itself (Fig. He). In this figure the loss of correla- 
tion between the electrons is evidenced by their escape 
in opposite directions. 
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V. TRIPLE AND HIGHER IONIZATION 



The model is easily extended to the case of simulta- 
neous removal of more than two electrons ||^. The key 
assumption is that the process is dominated by a sym- 
metric configuration of the electrons with respect to the 
field polarization axis. Specifically, we assume that all 
electrons move in a plane perpendicular to the field and 
that they obey a Cnv symmetry, which generalizes the 
C2v symmetry of the previous case. The reflection sym- 
metry limits the momenta to be parallel to the symmetry 
planes and thus confines the motion to a dynamically al- 
lowed subspace in the high-dimensional TV-body phase 
space. 

With the electric field directed along the z-axis the 
positions of the TV-electrons are Zi — z, pi — p and 
ifi = 27ri/iV, where (p^, (pi, z^) are cylindrical coordinates. 
The momenta of the electrons are all identical, pp^i — pp, 
Pz,i = Pz and Pip^i = 0. For this geometry the classi- 
cal Hamiltonian for N electrons, for zero total angular 
momentum along the field axis reads 

H{pp,p,,p,z,t) = N^^^^ + V{p,z,t), (4) 

with potential energy 

iV^ N(N-l) 
V = ~ . ^ + . \ . + NzFfit) cos{ut + c^) . 
^Jp^ + z^ 4psm(7r/7V) 

(5) 

The equipotential curves look very much like the ones 
shown for two particles and the dynamics is similar. One 
interesting aspect of this many electron configuration is 
that it is limited to at most 13 electrons: for larger num- 
bers of electrons the repulsion between overweights the 
attraction to the nucleus and no saddle configuration can 
be found. 



VI. FINAL REMARKS 

The present considerations suggest that correlated, 
non-seuqential multiple ionization in strong laser pulses 
proceeds through a saddle configuration with symmetri- 
cally moving electrons. The configurations can be seen 
analogous to the symmetrically escaping electrons in dou- 
ble ionization without field as discussed many years ago 
by Wannier. As in that case it is possible to derive a 
threshold law, which, however, is not only different from 
his but also much more difficult to observe because of the 
presence of the laser pulse. Further consequences of the 
model are under investigation. 
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